Radial quadrupole and scissors modes in trapped Fermi gases across the BCS phase 

transition 
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The excitation spectra of the radial quadrupole and scissors modes of ultracold Fermi gases in 
elongated traps are studied across the BCS superfluid-normal phase transition in the framework of 
a transport theory for quasiparticles. In the limit of zero temperature, this theory reproduces the 
results of superfluid hydrodynamics, while in the opposite limit, above the critical temperature, it 
reduces to the collisionless Vlasov equation. In the intermediate temperature range, the excitation 
spectra have two or three broad peaks, respectively, which are roughly situated at hydrodynamic and 
collisionless frequencies, and whose strength is shifted from the hydrodynamic to the collisionless 
modes with increasing temperature. By fitting the time dependent quadrupole deformation with 
a damped oscillation of a single frequency, we can understand the "jump" of the frequency of the 
radial quadrupole mode as a function of interaction strength which has recently been reported by 
the Innsbruck group. 

PACS numbers: 03.75.Kk,03.75.Ss,67.85.De,67.85.Lm 
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I. INTRODUCTION 

In the recent years, experiments with ultracold trapped 
fcrmionic atoms like 6 Li or 40 K have become a powerful 
tool for improving our understanding of strongly corre- 
lated Fermi systems. For example, the continuous tun- 
ability of the interatomic interaction via Feshbach res- 
onances made it possible to study the crossover from a 
Bose-Einstein condensate (BEC) of diatomic molecules 
to a BCS-type superfluid. The creation of vortex lat- 
tices in rotating Fermi gases with attractive interaction 
(scattering length a < 0) showed unambiguously that the 
superfluid BCS phase was indeed reached 

By studying the collective modes of these systems, one 
can obtain a lot of interesting information which cannot 
be obtained from static properties alone. Near the BCS 
phase transition, several collective modes of the trapped 
Fermi gases show interesting features like an extremely 
strong damping or a sudden change of the frequency 
@- H, HI H, Q- However, most of these experiments were 
done in the vicinity of the Feshbach resonance (/cfM ^ 1, 
where fop is the Fermi momentum in the center of the 
trap), and in this strongly interacting regime no theory is 
available so far which could quantitatively explain the ob- 
served damping effects. An obvious qualitative difference 
between the strongly and the weakly interacting regimes 
is that in the latter (fc^|a| < I) the normal phase is in the 
collisionless regime, whereas in the former (fci?|a| > 1) 
collisions can produce a hydrodynamic behavior of the 
normal phase [6[. 

The aim of the present paper is to discuss the proper- 
ties of experimentally accessible collective modes, namely 
of the radial quadrupole and scissors modes, in the BCS 
regime (fc^|a| < 1). In particular, we want to see how 
the properties of the modes change across the superfluid- 
normal phase transition. At zero temperature, the collec- 
tive modes can be described within the framework of su- 
perfluid hydrodynamics, but this is no longer true at non- 



zero temperature, where superfluid and normal compo- 
nents coexist. In Ref. [7[ it was proposed to use Landau's 
two-fluid hydrodynamics in this case. However, although 
this approach may be valid in the strongly interacting 
regime, it cannot be applied in the weakly interacting 
regime because of the lack of collisions. The theoretical 
framework which will be used in the present paper is the 
semiclassical transport theory developed for clean super- 
conductors by Betbeder-Matibet and Nozieres [8[ , which 
has recently been applied to the case of trapped Fermi 
gases 19, 101. In the case of a quadrupole oscillation in a 
spherical trap, this theory predicts a strong damping of 
the hydrodynamic mode in the superfluid phase due to 
its coupling to the collisionless normal component of the 
system, made of thermally excited quasiparticles. 

The article is organized as follows. Sec. HT1 gives a brief 
summary of the theoretical method (for details, see Refs. 
[§] and Qil). I n Sec. IIIII the results for the quadrupole 
mode are discussed, and it is shown that the frequcne; 
jump observed in a recent experiment at Innsbruck 
can be understood. In Sec. IIV1 results for the scissors 
mode are presented and Sec. [V] contains further discus- 
sions and conclusions. 



II. THEORETICAL METHOD 

The method used in the present work to describe col- 
lective modes in the BCS phase at finite temperature is 
based on the quasiparticle transport theory by Betbeder- 
Matibet and Nozieres [![. The basic degrees of freedom 
are the phase <j)(r,t) of the superfluid order parameter 
A(r, t), describing the collective flow of the Cooper pairs, 
and the quasiparticle distribution function v(r,p,t), de- 
scribing the normal- fluid component. In equilibrium, one 
can choose = 0, and the quasiparticle distribution 
function v reduces to the usual equilibrium Fermi func- 
tion f(E) = \/{e E / T + 1), where E = Vh 2 + A 2 with 
h = p 2 /2m + V(r) — /i, V and /i being the trap po- 
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tential and chemical potential, respectively. In the non- 



equilibrium case, one writes A 



, such that 



the collective velocity of the Cooper pairs is given by 
v s = —KV(f>/m. In this case, v is defined in the rest 
frame of the Cooper pairs, and the density p and current 
j are given by 



d 3 p f 1 h + 
(27rfi.) 3 V2 ~ 2El 
d 3 p 



{2nh) 3 m 



(1) 
(2) 



The subscripts + and — denote the time-even and 
time-odd parts of a given function (e.g., z/±(r,p, t) — 
(i/(r,p,t)±i/(r,-p,i))/2), 



n = \- V — nxp — n 



2m 



is the hamitonian and 



E 



h\ 



h- 



(3) 



(4) 



the quasiparticle energy in the rest frame of the Cooper 
pairs. 

The equation of motion of v can be written in the 
compact form 



(5) 



where {E, v\ denotes the Poisson bracket of E and v, 
similar to the usual Vlasov equation, which can be writ- 
ten as / = {h, /}, where / denotes the usual distribution 
function. The equation for <fi can be obtained from the 
continuity equation 



(6) 



which becomes a second-order differential equation for <J> 
if one inserts the explicit expressions for p and j . 

In order to describe collective modes, one writes V — 
Vq + Vi, v = vq + V\i and <f> = <pp + 4>i, with vp — 
f(E ) and (f)p = 0, the subscripts and 1 referring to 
equilibrium quantities and deviations from equilibrium, 
respectively. Assuming that the perturbation V± is weak, 
one keeps only terms which are linear in the deviations 
from equilibrium. The explicit equations of motion for </>i 
and v\ are quite cumbersome and can be found in Refs. 
@ and 0. 

Note that, contrary to Ref. the mean- field 

(Hartree) potential gp (g being the coupling constant) 
has been neglected in the definitions of h and h. If one 
wants to approach the crossover region (fcp|a| ~ 1), the 
expression gp for the mean-field potential is not applica- 
ble. One possible solution of this problem is to replace gp 
by the real part of the single-particle self-energy in lad- 
der approximation (Tlj . Anyway, previous experimental 
and theoretical studies showed that the mean-field shift 
does not modify the qualitative behaviour of the modes 



and leads only to minor changes of the mode frequen- 
cies. For example, in the recent Innsbruck experiment 
on the radial quadrupole mode the frequency shift 
in the collisionless normal phase from 2to r (u) r being the 
radial trap frequency) to ~ 2.luj r was attributed to the 
mean field. However, the mean field can considerably 
change the equilibrium density profile, thereby modify- 
ing the gap, the critical temperature, etc. For a quanti- 
tative study, the inclusion of the mean field is therefore 
desirable, but beyond the scope of the present paper. 

Except in very simple cases (e.g., a sound wave propa- 
gating in a uniform system Q), the equations of motion 
for v\ and 4>i have to be solved numerically. In Ref. [l(| , 
a test-particle method was proposed which is similar to 
the one often used for solving the Boltzmann equation. 
In the context of the linearized theory, it is useful to write 
V\ in the form 



df(E Q ) 

"i(r,P,t) = ,„ y\T,p,t) 
dEn 



(7) 



The advantage of this rewriting is that y is a smooth 
function in phase space, contrary to v\ [this is analogous 
to writing f\ = /o(l — fo)$ in the context of the lin- 
earized Boltzmann equation, which is a rather common 
trick] . The test-particle method consists in replacing the 
continuous function ^i(r, p, t) by a sum of S functions, 

ui{t, p, t) - yi (t) 6[r - K t (t)} S[p - Pi(t)] (8) 

i=i 

where yi(t) — 2/[Rj(t), Pi(t), t]. The trajectories of the 
test particles in phase space, (t) and P.; (t) , satisfy the 
following equations of motion: 



R 



dE 
dP 



dE 
<9R 



(9) 



These equations of motion do not conse rve the e nergy 
ho, but the quasiparticle energy Ep = y/jg + Aj). This 
can lead to surprising effects which make the numeri- 
cal solution more involved than that of usual Newtonian 
equations of motion. A well-known example is Andreev 
reflection: If a particle (h > 0) arrives at a point with 
A = E, it is reflected as a hole (h < 0), and vice versa. 

In the present calculations, the number of test particles 
is set to N v = 10 5 as in Ref. However, practically 

identical results are already obtained with N v = 10 4 . In 
fact, the large number of test particles was only needed 
in Ref. [l(| in order to obtain a stable mean field, while 
here, where the mean field is neglected, a much smaller 
number of test particles is sufficient. 

The time dependence of the coefficients y%{t) can rather 
easily be solved. For the time dependence of the phase, 
4> is expanded in a basis and the time dependence of the 
coefficients is determined by minimizing the violation of 
the continuity equation. In the case of the modes under 
consideration, the velocity field and therefore the form of 
(j) is well known. In this case, one or two basis functions 
are completely sufficient [lfj. 
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FIG. 1: Excitation spectra as function of the excitation fre- 
quency uj (in units of the radial trap frequency u r ) for different 
temperatures ranging from 0.157c to 0.75Tb (from bottom 
to top). For a better visibility of the strength contained in 
the narrow peaks at 2uj r , all spectra have been folded with a 
Lorentzian of width « 0.03tJ r . 



III. RADIAL QUADRUPOLE MODE 

Let us start by discussing the radial quadrupole mode. 
As it is usually the case in the present experiments, 
the unperturbed trap potential is supposed to be cigar- 
shaped: Vo(r) = |m[w r 2 (i 2 + y 2 ) + uo 2 z z 2 ] with u r > lo z . 
The radial quadrupole mode is excited at time t = by 
applying a perturbation of the form V±(r, t) = Vi(r)5(t) 
with Vi(r) oc x 2 — y 2 . The equations of motion are then 
solved for t > 0. The excitation spectrum is obtained 
as the Fourier transform of the time evolution of the 
quadrupole moment Q = (x 2 — y 2 ). 

We choose here parameters similar to the recent Inns- 
bruck experiments: N = 4 x 10 5 atoms, trap frequen- 
cies oj r = 2tt x 370 Hz, uj z = 2tt x 22 Hz. Since on 
the one hand our theory is only valid in the BCS limit, 
but on the other hand the scattering length of 6 Li sat- 
urates at quite a large value at high magnetic fields, we 
choose as a compromise l/k F a — —1.5. For these pa- 
rameters, the critical temperature within BCS theory is 
T c « 0.0587> « 43 nK (where T F = e F /k B is the Fermi 
temperature). However, one should remember that BCS 
theory cannot be trusted quantitatively. To give an ex- 
ample, the inclusion of the correlated density according 
to Nozieres and Schmitt-Rink keeping the number 
of atoms N = 4 x 10 5 and the parameter l/kpa = —1.5 
fixed, reduces the critical temperature by almost 20 % 
to T c ~ 0.049TV (k F and T F being defined as usual by 
the density of a non-interacting system with the same 
number of atoms). 

The resulting excitation spectra for different tempera- 
tures T are shown in Fig.[TJ At T — 0.15Tc, the spectrum 
is dominated by a peak situated at the hydrodynamic fre- 
quency, \2oj r . There is, however, a considerable broad- 
ening of the peak, i.e., damping, due to the coupling of 



the hydrodynamic mode to the thermally excited quasi- 
particles. At higher temperatures, the peak is further 
broadened and slightly shifted towards lower frequencies. 
At the same time, a second peak appears, which is situ- 
ated at 2io r , the frequency of the radial quadrupole mode 
in an ideal (collisionlcss normal- fluid) Fermi gas. Close to 
Tc, this mode is practically undamped (the width of the 
peak at 2u r in the uppermost spectrum in Fig. [T] is due 
to an artificial broadening which has been introduced in 
order to improve the visibility of the peaks). As one can 
clearly see, the peak position does not move from \[2io T 
to 2ui r with increasing temperature, but the strength is 
shifted from the lower to the upper peak. There seems 
to be no simple explanation for the downshift of the hy- 
drodynamic mode, but it can partially be understood as 
a consequence of the usual repulsion between two cou- 
pled modes. The general behavior is very similar to that 
found in previous work on the quadrupole mode in spher- 
ical traps 0, EE1 d ■ 

This finding suggests the following interpretation of the 
results of the Innsbruck experiment Q, where a sudden 
"jump" of the frequency as a function of the interaction 
strength was found. In the analysis of the experiment, 
the deformation of the atom cloud was fitted with an 
ansatz function containing only a single frequency and 
damping rate. Let us try what happens if we apply the 
same analysis to our theoretical results. 

In order to follow as closely as possible the experiment, 
we will do the calculation as a function of the coupling 
strength for fixed temperature, which we choose rather 
arbitrarily as T = 0.046TV. We also simulate the experi- 
mental initial conditions, which are different from the 5- 
like perturbation used above. In the experiment, the po- 
tential is first slowly deformed from Vo(r) to Vb(r)+Vi(r), 
with V\(y) rx x 2 — y 2 , giving enough time to the collisions 
to keep the momentum distribution spherical, and then 
the deformation Vi(r) is suddenly switched off at t = 0. 
Unfortunately there is one step in the experiment which 
cannot be simulated within the present linearized theory, 
namely the expansion of the cloud before the picture is 
taken. However, if one assumes that the picture taken 
after the expansion reflects in some way the state of the 
system at the moment when the trap was switched off, 
it seems plausible to assume that the expansion is not 
important for the observed frequency and damping, but 
only for the observed amplitude and phase of the oscilla- 
tion. 

In Fig. [5] the obtained time dependence of the 
quadrupole moment Q(t) is shown for different values 
of l/k F a. The solid lines are the calculations, while the 
dashed lines correspond to fits with functions of the form 

Q fit (t) = Ae- Kt cos(wt + </>) + Ce~ £t . (10) 

The parameters A, k, lo, 4>, C ', and £ are determined by 
minimizing the error 

((AQ) 2 ) = t~~ / dt {Q(t) - Q fit {t)f , (11) 

tmax JO 
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FIG. 2: Quadrupole moment as function of time for different 
values of 1/kpa ranging from —1.1 to —1.5 (from bottom to 
top) for fixed temperature T. 

where t max = 14 ms is chosen as in Ref. 5j. It can be 
seen that for l/fci?et = —1.1 and for 1/fcjra = —1.5 the 
fit reproduces almost exactly the full calculation, while 
for the intermediate couplings, in particular for 1/kpa = 
— 1.3, the fit is quite bad. 

In Fig. [3] the fitted frequency to and damping rate k 
are shown as functions of the parameter 1/kpa. The 
third panel shows the corresponding rms error of the fit, 
defined as ^ {(AQ) 2 ) /\Q(t = 0)|. There is a couple of 
striking similarities with the experimental results: At 
strong coupling, where Tq is large and hence T/Tc is 
small, the fit gives the hydrodynamic frequency. With 
decreasing coupling, i.e., increasing T/Tc, the damping 
increases and the frequency of the hydrodynamic mode is 
shifted downwards. The quality of the fit becomes quite 
bad (this corresponds probably to the shaded area in Ref. 
(H), and at some critical coupling strength the fit jumps 
to the collisionless (normal-fluid) mode. As seen in the 
experiment, the damping rate of the collisionless mode is 
weaker than that of the hydrodynamic mode before the 
jump. 

It should be noticed that with the present choice of 
parameters, the jump happens at 1/kpa ~ —1.32, cor- 
responding to T/Tc ~ 0.6. This means that in the cen- 
ter of the trap the system is still superfluid, but since 
this region contributes only very little to the quadrupole 
moment, the quadrupole mode behaves almost as if the 
system was in the collisionless normal-fluid phase. 

In spite of the surprising qualitative agreement with 
the findings of Ref. [f|, there is no exact agreement 
on a quantitative level. First of all, in the experiment 
the sudden transition from hydrodynamic to collisionless 
normal-fluid behavior happens already at a stronger cou- 
pling than in the calculation. The position of the jump 
depends on the temperature, and if the temperature in 
the experiment had been higher than ours (0.046Tf), this 
could explain the difference. In addition, one should re- 
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FIG. 3: From top to bottom: Fitted frequency u>, damping 
rate k, and corresponding rms error of the fit (definition see 
text) for the radial quadrupole mode as functions of 1/kpa 
for fixed temperature T. The dotted lines in the upper panel 
indicate the prediction of superfluid hydrodynamics (\ /f 2uj r ) 
and the prediction for an ideal Fermi gas (2ui r ). 

member that the critical temperature predicted by BCS 
theory is too high. The second problem is the damping 
rate: The experimental damping rates are larger than the 
theoretical ones by roughly a factor of 2 near the jump, 
and even worse in the case of the collisionless normal- 
fluid mode far away from the transition. In fact, the 
experiment is not in the collisionless limit and collisions 
still play a role. In order to improve the theory in this 
respect, one would have to include a collision term sim- 
ilar to that of the Boltzmann equation into the quasi- 
particle transport equation. The small anharmonicity of 
the experimental trap potential could be an additional 
source of damping. However, a calculation where the ra- 
dial harmonic potential, ^muj 2 (x 2 + y 2 ), was replaced by 

a Gaussian one, — Uoe~( x +y ^ d , with a realistic width 
d = 0.05 mm, showed that this effect can explain only 
~ 10% of the experimentally observed damping in the 
normal phase. The third point concerns the frequencies: 
On the one hand, although the theory gives a down-shift 
of the hydrodynamic mode before the transition, it is not 
as strong as the down-shift observed in the experiment. 
On the other hand, the theoretical frequency of the col- 
lisionless normal-fluid mode above the transition is too 
low, because the mean-field has been neglected. 



IV. SCISSORS MODE 

In the case of the radial quadrupole mode the velocity 
field is always irrotational and it does not change signifi- 
cantly from the superfluid to the normal-fluid phase. The 
different frequencies are a consequence of the Fermi sur- 
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FIG. 4: Excitation spectra of the scissors mode as functions of 
the excitation frequency uj for different temperatures ranging 
from 0.3Tb to 0.9Tb (from bottom to top). The spectra have 
been folded with a Lorentzian of width ~ 2n x 10 Hz. The 
arrows indicate the predicted frequencies in the hydrodynamic 
(v/a>! + My) and collisionless limit (u) x ±ui y ), respectively 



face distortion during the oscillation in the collisionless 
normal-fluid phase. This is different from the "scissors 
mode," which has also been studied in a recent exper- 
iment [6J. In this experiment, the trap was deformed 
in the xy plane: Vb(r) = ^m[uj 2 x 2 + Lo^y 1 + oj 2 z 2 ], 
with lo x > Wy S> w z . The scissors mode corresponds 
to a motion where the deformation of the atom cloud 
in the xy plane is not aligned with that of the trap, so 
that the angle between them oscillates. In the normal- 
fluid phase, such a motion can be realized by rotating 
the gas back and forth around the z axis [velocity field 
v(r) oc e z X r], but also by an irrotational motion [veloc- 
ity field v(r) oc V(xy)]. In an ideal gas, these two kinds 
of motion correspond to two distinct modes with frequen- 
cies iO- and u>+ , respectively, given by tu± = uj x ±uj y . In 
the superfluid phase, of course, the rotational mode is 
suppressed and only the irrotational one exists. In the 
zero-temperature limit, its frequency approaches the hy- 
drodynamic value uj sc = ^Jui 2 . + ujy [ll|. 

The excitation of the scissors mode consists in turning 
the trap potential around the z axis by a small angle. 
In the small-amplitude limit, this corresponds to an ex- 
citation operator Vi(r) oc xy. The excitation spectra of 
the scissors mode for different temperatures, which are 
obtained analogously to those of the radial quadrupole 
mode (i.e. Fourier transform of (xy) after a J-like per- 
turbation), are shown in Fig.0]for the case of N = 3x 10 5 
atoms in a trap with oj x — 2it x 580 Hz, tv y — 2n x 270 
Hz, uj z = 2tt x 22 Hz, and l/k F a = -1.5. At T = 0.3T C , 
the spectrum has a single broad peak at the hydrody- 
namic frequency. The damping rate (width) is compa- 
rable with that of the hydrodynamic radial quadrupole 
mode at T = O.STq. With increasing temperature, the 
two collisionless normal-fluid modes appear at the pre- 
dicted frequencies, while the strength of the hydrody- 



namic mode decreases and its frequency is slightly shifted 
downwards. At T = 0.6Tc and 0.75Tc the scissors-mode 
spectrum has three pronounced peaks. There is a strik- 
ing analogy between these results and the experimental 
observations concerning the scissors mode in a BEC with 
a thermal cloud, where the BEC oscillates with uj sc which 
is slightly shifted downwards while the thermal cloud os- 
cillates with uj± [161 ]. 

The results of Fig. 0] are quite different from those 
found in the experiment [6J. However, this is not sur- 
prising, since the experiment was not done in the BCS 
limit but for 1/kpa = —0.45, where in particular the col- 
lisions are very important. As one can see in Fig. 2 of Ref. 
[(|, the phase transition in the experiment does not lead 
directly from the superfluid to the collisionless regime, 
but the system enters first a collisionally hydrodynamic 
regime (i.e., the mode frequency does not change) and 
reaches the collisionless phase only at much higher tem- 
perature. As shown in Ref. (l7| . the transition from colli- 
sional hydrodynamic to the collisionless regime leads to a 
continuous increase of the upper (irrotational) branch of 
the scissors mode. This explains why in the experiment 
no sign of a sudden jump of the scissors mode frequency 
has been observed. It would be interesting to see if in the 
BCS regime, e.g. at 1/kpa = —1.5, the experiment would 
observe a similar jump as in the case of the quadrupole 
mode. 



V. CONCLUSIONS 

In the present paper, the behavior of the radial 
quadrupole and of the scissors mode across the BCS 
phase transition was theoretically studied in the frame- 
work of a semiclassical theory, in which superfluid hydro- 
dynamics is coupled to a quasiparticle transport theory 
for the normal component. The theory proved to be ap- 
plicable to realistic situations (large numbers of atoms, 
arbitrary trap geometry), and it provides an explana- 
tion of the surprising behavior of the radial quadrupole 
mode as a function of 1/fcpa as observed by the Innsbruck 
group [5]: The damping near the BCS phase transition is 
due to the coupling between the hydrodynamic mode and 
the thermally excited quasiparticles. In addition to the 
damping, this coupling leads to a two-peak spectrum, the 
two peaks belonging, respectively, to the hydrodynamic 
mode of the superfluid component and the collisionless 
mode of the normal-fluid component. The sudden jump 
from the hydrodynamic to the collisionless frequency re- 
ported in Ref. _5J is probably a consequence of the fit 
with a single frequency in the analysis of the experiment. 

It should be mentioned that an alternative explana- 
tion of the damping was suggested in Ref. [18j: The 
frequency shift and damping are not necessarily a sig- 
nal for a superfiuid-normal phase transition. They could 
also be related to the breakdown of superfluid hydro- 
dynamics when fko is of the same order of magnitude 
as the pairing gap A. In Ref. effects of this kind 
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were quantitatively studied by using a fully quantum me- 
chanical formalism [quasiparticle random phase approx- 
imation, (QRPA), corresponding to the time-dependent 
Bogoliubov-de Gennes equations in the case of small de- 
viations from equilibrium], and it was shown that in- 
deed, even at zero temperature, considerable damping 
and deviations from the hydrodynamic frequency can oc- 
cur if fajjj A is not small. Although such effects may also 
play a role in the recent experiments (maybe they can 
give rise to an additional downshift of the hydrodynamic 
quadrupole mode, which would help to quantitatively ex- 
plain the data) , they are probably less important than the 
coupling between the hydrodynamic mode and the ther- 
mally excited quasiparticles considered here: In Ref. 
it was shown that for a ratio of mode frequency to zero- 
temperature gap of hw I A(T = 0) ~ 0.3 (which is similar 
to the most unfavorable ratio of all the cases considered in 
the present paper), the results of the semiclassical the- 
ory and of the quantum mechanical QRPA calculation 
for the quadrupole mode in a spherical trap are in good 
agreement, even if the finite-temperature gap is smaller. 

Possible extensions of this theory which are needed for 
giving more quantitative results at the limit between the 



BCS and the crossover regime (fcpH ~ 1) have been 
mentioned. These are, on the one hand, the inclusion 
of quasi-particle collisions in order to reproduce the ob- 
served damping, and on the other hand the inclusion of 
the single-particle self-energy in order to account for the 
"mean field" energy shift. 

However, these extensions may not be enough to give a 
complete description of the BCS-BEC crossover (/cfM <L 
1), including the unitary limit (/cfM - ► oo). In this 
regime, the normal-fluid component of the system is not 
only made of thermally excited quasiparticles, which are 
included in the present theory, but also of pair fluctua- 
tions [l!|, which persist even above the critical temper- 
ature Tc up to the "pair breaking temperature" T* and 
which cannot easily be accounted for in the present the- 
oretical framework. 
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